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Abstract 

We  consider  the  problem  of  modeling  dynamical  effetcs  of  impact  of  an 
elastic  body  on  a  flexible  beam.  We  derive  a  nonlinear  integral  equation  by 
using  the  Hertz  law  of  impact  in  conjunction  with  the  beam  equation.  This 
equation  does  not  admit  a  closed  form  solution.  We  demonstrate  the  existence 
of  solutions,  derive  a  reliable  numerical  method  for  computing  solutions,  and 
compare  the  numerical  results  with  those  obtained  by  others. 


1  Introduction 

High  precision  control  of  robotic  manipulators  has  been  becoming  increasingly  im¬ 
portant  in  a  variety  of  applications,  e.g.  laser  beam  technology,  semiconductor  safer 
manufacturing  etc.  This  requires  paying  extra  attention  to  the  usual  dynamical 
effects  as  well  as  taking  into  consideration  otherwise  ignored  features  such  as  dy¬ 
namical  effects  due  to  impact.  This  paper  focuses  on  the  latter  aspect. 

For  the  sake  of  simplicity,  we  only  consider  an  elastic,  beam  subject  to  impact 
forces  occuring  from  contact  with  an  elastic  body.  Here  we  restrict  attention  to  the 
problem  of  modeling,  existence  of  solutions  to  the  model,  and  the  computational 
aspects.  Issues  such  as  how  to  control  the  manipulators  to  minimize  the  spurious 
effetcs  due  to  impact  will  be  addressed  in  the  future. 

Numerous  attempts  have  been  made  to  accurately  model  dynamical  effects  of 
impact  in  robotics  oriented  applications  in  the  recent  years[l,2,3j.  Consideration  of 
displacement  and  use  of  Hertz’s  law  of  impact  at  the  region  of  contact  seems  to  be 
the  most  successful  a.pproach[4] .  When  the  contact  involves  a  flexible  beam,  Hertz’s 
law  of  impact  leads  to  a  nonlinear  integral  equation  called  the  Hertz  equation,  which 
incorporates  the  effects  of  local  elastic  deformation  at  the  region  of  contact [5].  This 
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model  has  been  widely  applied  to  various  impact  situations,  and  the  experimental 
results  obtained  in  [2] [6]  well  support  the  validation  of  this  equation.  Unforuntely, 
this  nonlinear  equation  does  not  admit  a  closed  form  solution.  Timoshenko[4]  used 
small-increment  method  to  give  numerical  solutions,  and  it  became  the  basis  for 
evaluting  other  approximation  methods.  Some  other  approximation  methods  also 
give  very  satisfactory  results.  One  of  them  is  the  application  of  the  energy  method 
devised  by  Zener  and  Feshbach[8],  and  applied  by  Lee[5]  to  central  impact  of  a 
sphere  on  a  simply  supported  beam.  These  approximation  methods  were  proposed 
without  establishing  the  convergence  or  even  existence  of  a  unique  solution. 

In  the  section  2  of  this  paper  the  impact  problem  is  properly  formulated  and  the 
Hertz  equation  is  derived  through  Hertz  law  of  impact.  In  section  3  we  will  discuss 
some  basic  properties  of  the  Green’s  function  associated  with  the  Euler  Bernoulli 
beam  equation.  This  equation  is  used  to  model  the  motion  of  the  beam.  In  section 
4,  we  will  establish  the  existence  and  uniquness  of  solutions  of  the  Hertz  equation 
by  applying  the  contraction  mapping  principle.  In  section  5,  a  numerical  technique 
based  on  the  contraction  mapping  principle  is  presented.  Various  examples  are 
discussed,  and  the  results  obtained  by  different  approaches  are  compared. 

2  Formulation  of  the  Problem 

For  our  purpose,  the  impact  problem  can  be  formulated  as  follows;  a  beam  is  struck 
transversely  by  a  mass  m  having  a  spherical  surface  at  the  point  of  contact  with 
initial  moving  velocity  vq.  We  further  assume  that  the  Hertz  law  of  impact  is  valid 
,  i.e, 

«  =  A-[/(*)]2/3,  (2.1) 

where  a  is  the  relative  approach  of  the  striking  body,  f(t)  is  the  contact  force  and  K 
is  the  Hertz  c.onstant[7],  which  is  determined  by  the  local  geometry  (i.e.  the  radius 
of  curvature)  and  material  properties  of  impacting  objects.  Relative  approach  is  the 
difference  between  the  displacement  of  the  beam  and  the  contacting  body,  measured 
from  instant  of  initial  contact.  Hence, 

a  =  s(t)  -  w(x *,t),  (2.2) 

where  w(x*,t )  is  deflection  of  the  beam  at  the  point  of  contact  x*,  s(t)  is  the 
displacement  of  the  ball  under  of  the  contact  force  f(t),  and  is  given  by, 

•s(t)  =  vQt  -  —  f  f(r)(t  -  T)(lr.  (2.3) 

in  Jo 

From  the  equations  (2.3)  ,  (2.1)  and  (2.2),  we  obtain  the  nonlinear  integral  equation, 

K[f(t)]2/3  =  Vo t  ~  ~  f  -  r)r/r  —  w(x*,t).  (2.4) 

in  Jo 
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Figure  1:  Sketch  of  the  Displacement 

Due  to  the  nonlinear  term  A'[/(<)]2/3,  it  is  impossible  to  obtain  a  closed  form  solution 
for  (2.4).  Before  we  can  carry  out  further  analysis,  it  is  necessary  to  represent  the 
deflection  of  the  beam  at  the  point  of  contact. 

3  Deflection  of  the  Beam  and  Green’s  Function 

If  we  only  restrict  to  transverse  vibrations  and  assume  that,  the  beam  is  long  and 
slender,  the  transverse  shear  and  torsional  effects  may  be  neglected,  and  the  dynam¬ 
ics  of  the  beam  can  be  described  by  Euler  Bernoulli  beam  equation.  The  deflection 
of  the  beam  is  in  turn  obtained  by  solving  the  following  partial  differential  equation, 

+  El  A  w  =  f(x,t)  0  <  x  <  /,  (3.1) 

dir 

in  which  A  =  l  is  the  length  of  the  beam,  p  is  the  mass  density  and  El  is  the 
bending  stiffness  (here  p  and  El  are  assumed  constants),  /(:»:,  t)  is  the  distributed 
load.  The  deflection  of  the  beam  is  uniquely  determined  when  equation(3.1)  is  solved 
under  the  appropriate  initial  and  boundary  conditions.  If  we  assume  that  the  beam 
is  at  rest  just  before  impact,  the  initial  conditions  are, 

t/>(:r,0)  =  w(x,  0)  =  0.  (3.2) 

Various  boundary  conditions  of  interest  can  be  described  as 

B;w(x,t )  =  0,  i  =  1,2,  x  =  0,/.  (3.3) 

where  B-,  is  a  linear  homogeneous  differential  operator  of  maximum  order  3. 
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The  concentrated  load  /(f)  is  obtained  as  a  limiting  process  of  a  uniformly 
distributed  load  f(x ,  t)  over  a  small  range  26  of  the  beam.  Thus,  by  letting  f(x,  t )  — ► 
oo  while  6  — >  0,  the  contact  force  /(f)  is  obtained  by, 

rx*+S  _ 

/(f)  =  lim  /  f(x,i)dx  (3.4) 

6 — ►(),/  — *oo  Jx*—6 

One  of  the  most  popular  methods  for  analyzing  partial  differential  equations  is 
the  Green’s  function  method[9].  With  the  aid  of  a  Green’s  function,  the  solution  of 
a  certain  class  of  PDE  can  be  expressed  as  an  integral  . 

Definition  3.1  A  function  G(x,(;t)  £  L2[0,l]  is  called  a  Green’s  function  of  (3.1)- 
(3.3),  if  it  satisfies  the  following  conditions: 

i)  vis  a  function  of  the  argument  x,  it  satisfies  the  homogenous  differential  equa¬ 
tion  ,  i.c.  f  =  0,  everywhere  except  at  x  =  (  where  it  may  have  a  singularity. 

ii)  As  a  function  of  the  argument  t,  it  satisfies  the  homogenous  differential  equa¬ 
tion  everywhere  except  at  t  =  0  where  it  may  have  a  singularity. 

Hi)  .4 .s  a  function  of  the  argument  x  ,  G(x,(;t)  satisfies  the  boundary  condition 
(3-3). 

iv)  For  the  initial  condition  ,  it  satisfies 

G(x,(;  0+)  =  0,  and  =  (3.5) 

Note  that  (i)  and  (ii)  above  leads  to, 

{Ei  A  =  HoOHt).  Ot.fi) 

Since  the  PDE  (3.1  )-(3.3)is  self-adjoint,  G{x,  (;  t)  is  symmetric  with  respect  to 
x  and  (,  and  can  be  expressed  in  terms  of  eigenfunction  expansition.  Therefore,  it 
can  be  shown  that  the  Green’s  function  for  the  PDE(3.1)-(3.3)  can  be  expressed  in 
the  form  (we  refer  the  reader  to  [10]  for  details), 

G(x,G,t)  =  f^Wk(x)Wk(O^^H(tfi  (3.7) 

Li  Wk 

where  H(t)  is  the  unit  step  function,  {Wk{x)}rlf=l  is  an  orthnormal  basis  of  eigen¬ 
functions  and  j  are  the  corresponding  eigenvalues.  It  is  easy  to  show  that  a 

representation  of  the  solution  to  the  PDE  (3.1)-(3.3)  in  terms  of  the  Green’s  function 
G(x,(-,t)  is 

w(x,t)  =  /  /  G(x,C;t- T)f(C,T)d(d.T.  (3.8) 

Jo  Jo 
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For  the  impact  problem,  since  the  contact  can  be  treated  as  point  contact,  con¬ 
tact  force  has  the  special  form  (3.4).  Hence  equation  (3.8)can  be  further  simplified 
as 

w(x,t )  =  f  G(x,x*;t  -  T)f(r)dT.  (3.9) 

Jo 

For  simplicity,  we  write  G(x*\t)  instead  of  G(x*,x*;t)  in  the  rest  of  the  paper. 
From  the  equations  (2.4)  and  (3.9),  the  Hertz  equation  will  be 

A'[/(0]2/3  =  v0t-—[  f{r)(t  —  t)(Jt  —  f  G(x*;t  -  r)/(r)dr.  (3.10) 

Til  Jo  Jo 

4  Analysis  of  the  Hertz  Equation 

Though  the  Hertz  equation  (3.10)  hasn’t  been  analyzed  in  any  great  detail  in  the 
existing  literature,  some  approximation  methods  for  solving  it  have  been  presented 
in  some  detail.  Our  viewpoint  is  that  theoretical  analysis  is  necessary  for  both 
proving  the  validation  of  this  equation  and  developing  efficient  numerical  methods. 
Contraction  mapping  technique  is  employed  here  to  show  that  a  unique  solution 
exists  for  the  Hertz  equation.  Before  invoking  the  contraction  mapping  theorem, 
some  simplifications  are  necessary.  Let, 

L(t)  =  t  +  mG{x*]t).  VI  >  0  (4.1) 

Equation  (3.10)  can  be  rewritten  as, 

/(*)2/3  =  %t  -  ~  f  f(r)L(t-r)dr,  (4.2) 

in  j  o 

where  v'0  =  cq / K ;  m!  =  mK; 

f(t)  =  [v'0t  -  —  (  f{r )L(t  -  t)(J.t]3/2  (4.3) 

=  K*  -  ~1  I  /(r)(^  -  T)dT  -  T  f  /(rKT'(**5 1  ~  T)dT\i/2-  (4-4) 

m  Jo  Jv  ./o 

Note  that  is  assumed  to  be  positive  always.  Both  equations  (4.3)and  (4.4)will  be 
used  in  the  following  analysis. 

Theorem  4.1  (Contraction  Mapping  Theorem)  Let  X  be  a  Banach  space,  and  B  be 
a  closed  subset  of  X.  Let  P:  B  — >  B  be  an  operator  satisfying  the  following  condition: 
3  p  <  1  such  that 


\\Px-Py\\  <  p\\x  —  j/||,  V.t ,  y  G  B. 
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Then 

a)  P  has  exactly  one  fixed  point  in  B  (denoted  by  x*). 

b)  For  any  x0  G  B,  the  sequence  { xn }°°  defined  by 

®n+l  —  P^rii  n  >  0 

converges  to  x*.  Moreover, 

ll*n-*1l  <  t^-HPxo-xoII. 

1-  p 

A  proof  of  this  well  known  theorem  can  he  founded  in[llj.  We  will  use  this  theorem 
as  the  main  tool  to  show  that  equation  (4.3)  has  a  unique  solution  by  constructing 
a  contraction  operator  P  on  an  appropriate  closed  subset  B  of  a,  Banach  space. 

Theorem  4.2  Suppose  that  the  Green’s  function  G(x*\t)  is  uniformly  bounded,  over 
[0,/].  Then  there  exists  a  small  enough  S  >  0  such  that  (f.3)  has  a  unique  continu¬ 
ous  solution  for  t  £  [0,  A], 

Proof  :  Let  M  >  0  be  such  that, 

\G(x*;t)\  <  M  Vi  >  0  and  Vx*  €  [0,/j; 

Let  NX)  be  a  sufficiently  large  constant.  Let  A  >  0  be  small  enough  such  that, 

W)  ^  ~  [l/2m'+M/2K]  ’ 

(Hi)  2(A2/m'  +  M6/K)^(v'a  +  MN/I<)S  <  1. 

Our  Banach  space  here  is  C  [0,  A],  the  space  of  continous  real  valued  factions  from 
[O’  A]>  endowed  with  the  sup  norm,  i.e.  \\f IU  =  sup<ef0  « \f(t)\.  Let  us  define  the 
mapping  P  :  C'[0,A]  C[0,A]  by, 

1  P  3/2 

p/(0  =  W0t  -  —  /  f(T)L(t  -  r)dr]  ,  Vte[0,A].  (4.5) 

The  domain  of  P  is  defined  by, 

B[0J}=  {/(•)€  C[0,  AJ;  A  > /(•)>(); 

1 

va 1  ~  ~ 7  /  f(T)W  -  r)dr>  0  V*€[0,A]}.  (4.6) 

lit  J o  ' 

Obviously,  /?[(),  A]  is  a  closed  subset  of  the  Banach  space  of  continuous  functions  on 
[0,A], 
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The  rest  of  the  proof  is  divided  into  two  parts:  first,  we  show  that  P  maps  i?[0,  <*>] 
into  itself.  Then  we  show  that  P  is  a  contraction  mapping  on 

a)  Pf  >0  V/  €  B[0,  <$]  by  definition.  Let  us  define  mapping  F  :  B[0,  <5]  — > 
B[0,S]by, 

Ff(t)  =  ivot  !  f(T)L(t  -  r)dr ] 
m'  Jo 

f  g  £[0,tf]=>-  £?fof(r)(t  -  t )dr  >  0.  Hence, 

\Ff(t)\  =  I  Wot - l-  f  f(r)(t  -  r)dr  -  ^  f  f(r)G(x*]t  -  r)dr]\ 

nr  Jo  li  Jo 

<  v’o t  +  jr  j  \f(T)\\G(x*;t  —  r)\dr 

<  (v'o  +  MN/K)t. 

Ff(t )  >  0  by  definition.  Therefore, 

Ff(t )  <  (v'o  +  MN/K)t=>Pf(t)  <  [K  +  MN/K)t]3'2  <  t,  Vi  €  [ (M]. 

FP.f(t)  =  v'0 1  -  —  f  Pf(r)(t-r)dT  -  ~  [  Pf(r)G(x*;t  -  r)dr 

ni  Jo  K  Jo 

>  vo t  -  — t  [  (t  -  T)rdr  -  77  f  tM dr 

m  Jo  A  Jo 

.  i  t2  t2M 
-  Vl)t  ~  m>  2  "  1 Y 
>0  Vi  G  [0, i>] 

Thus,  we  have  shown  that  PB[0,8]  C  B[0,  i>]. 

b)  Vfi,f2eB[ 0,  A], 

Pfi(t)  -  P.h(t)  =  Wot  -  —  f  /i(r )L(t  -  r)r/r]3/2  -  [«£t  -  —  I  /2(r)Z(i  -  r)drf/2 

m  Jo  m'  Jo 

Let  x(t)  =  [v'ot.  -  —  f  fi(r)L(t-  r)4r]1/2 

nv  Jo 

y(t)  =  Wot-^jf  /2(r)G(i  -  r^r]1/2 
m 

5mce  /i,/2  6  J9[0,tf],  =>  .r(i)  >  0,  v/(i)  >  0  Vi  6  [0, tf] 

A/i(i)  -  Pfi(t)  =  x3(t)  -  ?/(i)  Vi  €  [0, S] 

\P.fi(t)  -  Ffi(t) |  <  |.r2(i)  -  i/2(i)||.r(i)  +  sf(i)| 
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*2(0  -  y2(t) 


\x\t)  -  y\t) | 


=  ~  ~i  J0  “  r)rfr  -  ('’(V  ~  ~i  JQ  Mr)L(t  ~  T)dT ) 

=  ^7  Jq{Ht)  ~  hix))L{t  -  T)d.r 

=  l~7  ih(T)  ~  fi(T))L(t  —  r)dr\ 

<  ~~i  /  |/a(r)-/i(r)|(t-r)dr+-^  /  |/a(r)  -  /i(r)||G(x*;  *  -  r)|rfr 

m'  Jo  A  Jo 

<  (*2/m'  +  Mt/JOII/a-/i||oo 

<  (AVm' +  MA/ A)!  I/2-/1 1 U 


K0  +  y(0l  <  KOI  +  M*)l  <  2^0*  +  MJV///1  <  2^/(«'  +  Mtf/A' )$ 
|P/i  W  -  P.m  |  <  |*2(0  -  y2(<)IN*)  +  2/(01 

<  2(62/m'  +  MA/A)^(u' +M7V//OA||/2(-)  -  /i(-)IU 

<  p\\f2  ~  /llloo  V2  G  [0,  ^)] 

where  p  =  2(S2/m'-\-Mf>/K)yJ(v'Q  +  MN/K)f>,  and  by  the  property  (in)  of  A,  p  <  1. 
Therefore, 


\\Pf1-Pf2W00  =  sup  |P/i(t)  -  P/2(i)| 

<e[o,«] 

<  HI/2  -  /llloo 

so  that  P  is  a  contraction  mapping  on  I?[0,A]. 

Finally,  using  the  theorem  /,.J.  it  follows  that  the  mapping  P  has  a  unique  fixed 
point  in  /?[(),  A] .  It  is  clear  that  f  is  a  solution  of  the  equation (Jt.  3)  over  [0,A]  iff 
P f  =  f,  i.e.  f  is  a  fixed  point  of  P  otter  /J[0,  A].  This  completes  our  proof. 

The  above  theorem  shows  that,  a  unique  solution  exists  over  t  6  [0,  A]  for  some 
small  6.  Our  interest  is  to  find  the  impact  force  variation  during  the  entire  contact 
period.  The  following  theorem  will  establish  this  global  result. 

Theorem  4.3  Suppose  that  the  equation  (f.3)  has  a  local  unique  solution  over  [0,  A] 
for  some  sufficiently  small  A.  If  /(A)  >  0,  then  3e  >  0,  such  that  the  equation (Jf. 3) 
has  a  unique  solution  on  [0,  A  +  e]. 

Proof  :  The  argument  is  similar  to  the  proof  of  the  local  version.  We  will  only 
carry  out  details  of  a  crucial  step  here. 
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Let  g  :  [0,£]  — >  R  be  the  unique  solution  of  (4-3)  established  in  the  proof  of  the 
theorem  4-3.  Let  N  be  a  positive  number  larger  than  ||s||oo-  Let  e  >  0  be  a  small 
positive  constant.  Let 


■D[0, 6  -f-  f]  — 


{/  £  C[Q,6  +  ej;  /  |[0^] 

v'0t  -  — -  f  /(r )L(t  -  t )dr  >  0 
m  Jo 


Vf  €  [0, 6  +  f]}- 


(4.7) 


Clearly  D  is  a  closed  subset  of  (C[0,8  -f  e],  ||  •  ||oo)- 


Let  F  :  jD[0,  6  +  e\  — >  C[0, 8  +  e]  be 

Ff(t)  =  v'0t  -  ^7  /  f(r)L(t  -  r)d,T, 
nv  Jo 

and,  let  P  =  2 . 


We  will  show  that  for  small  enough  e,  P(D)  C  D,  and,  P  is  a  contraction 
mapping,  thus  establishing  the  theorem.  Note  that  it  follows  easily  as  in  the  proof  of 
theorem  4-2  that 

\Pf(t)\  <  N  \/te[0,S  +  e}. 

The  crucial  step  is  to  show  that  FPf(t)  >0  V/  £  [0,  8  +  e],  V/  £  D.  Now, 


FPf(t)  =  /(/)  \/t  <  8  since  f  |[0i^  satisfies  (4-3).  For  0  <  t  <  e, 

FPf{t  +  8)  =  v'08  -  —  I  Pf(r)(t  +  8  -  r)d,T  -  77  /  Pf(T)G(x*;  t  +  8  -  r)d,T 

in'  Jo  A  Jo 

+v'ot  ~  ~1  I  Pf(T  +  b)L(t  ~  T)dr 
m  Jo 

=  Vu 8  -  -7  /  Pf{r)(8  -  t )dr  -  ^7  /  Pf(r)G(x *;  8  -  r)d.T 
m  Jo  A  Jo 

I  f{T )dr  -  77  /  f(r)(G(x*,t  +  8  -  r)  -  G[x*\ 8  -  r))dr 
111  Jo  A  Jo 

+'"o<  -—if  Pf(T  +  b)L(t  -  r)d,T 
nv  Jo 

=  [f2/\b)-—  f*  f(T)dr-G(x*;t)  +  G(x*;0) 

rn  Jo 

pigt  -  —  f  Pf(r  +  8)L(t  -  T)dr] 
in  Jo 

where,  I)  =  Jo  .f(T)G(x*G-  +  8  —  r)dr.  By  the  continuity  of  the  Green's 

function,  for  small  enough  e, 

\G(x*,t)  -  G{x* ,Q)\  <  l-pl\8)-  VfG[ 0,*] 
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For  such  e 


FPf(t  +  S)  >  0  V/€[0,e]. 

Hence,  we  have  shown  that  PD[ 0,  S  +  e]  C  .D[0,  S  +  e]. 

It  is  easy  to  show  that  P  is  a  contraction  mapping.  Existence  of  a  unique  solution 
on  [0, 6  +  e]  follows  at  once  now. 

Remark  4.1  The  global  version  has  the  following  physical  interpretation.  The  con¬ 
dition  f(t )  >  0  means  that  the  contact  is  in  progress  at  t  >  0;  finally,  f(T)  =  0 
means  that  the  objects  are  just  about  to  cease  to  be  in  contact,  i.  e.  T  is  the  impact 
duration. 


5  Numerical  Methods 


Some  approximation  methods  to  solve  the  equation  (4.3)  have  been  proposed  in  the 
literature.  One  of  them  is  the  energy  method[5].  It  is  simple  and  fast,  although 
less  accurate  than  the  small-increment  method.  Hence,  the  results  obtained  by 
the  energy  method  can  be  used  as  a  good  initial  condition  for  the  more  accurate 
methods.  The  gerenal  solution  by  the  energy  method  has  a  sinusoidal  form: 

fo{t)  =  Kosin(xt/T0)  V/  €  [0,To], 


where,  T0  is  the  impact  duration  calculated  from  the  energy  method,  and  K 0  is  a 
constant.  Inspired  by  the  contraction  mapping  theorem,  we  develop  a  numerical 

method  using  the  successive  Picard  approximations;  /o,  Pfo,  PP.fo , . ,  where  the 

initial  condition  fo  is  obtained  from  the  energy  method,  and  P  is  the  contraction 
operator  defined  by, 


P.t\l)  =  K*  -  77  /  f(r)L{t  -  r)dr] 

nr  j o 

Closed  form  solution  of  the  first  iteration  is, 

fi(t)  =  K< - ~7  /  fo{r)L(t  -  t )dr] 

nr  Jo 

Numerical  results  are  given  below  for  two  examples. 


3/2 

? 


3/2 


v/e  [o,r0]. 


v*  e  [o, To]. 


In  example  i,  impact  occurs  at  the  center  of  a  simply  supported  beam.  Fig.  1 
shows  the  solution  obtained  from  the  energy  method,  and  gives  the  force  magtitude 
and  duration  very  close  to  that  obtained  by  using  small-increment,  method,  although 
the  force  shape  differed  in  some  ways.  In  order  to  apply  the  Picard  approximation 
method  to  this  case,  we  need  first,  to  show  that  the  Green’s  function  associated  with 
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Force(Newton) 


Figure  2:  Central  impact  on  a  simply  supported  beam 


Figure  3:  Tip  impact  on  a  cantilevered  beam 
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this  case  is  uniformly  bounded.  Details  are  given  in  Appendix  A.  The  result  obtained 
by  3-step  picard  approximation  is  very  good,  and  the  computational  complexity  is 
just  1/3  of  the  small-increment  method. 

In  example  2,  impact  is  at  the  tip  of  a  cantilevered  beam.  We  see  that  fairly  large 
errors  occured  by  using  the  energy  method.  On  the  other  hand,  the  new  numerical 
method  gives  excellent  results  even  after  just  one  iteration,  which  has  a  closed  form 
solution.  Fig. 3  also  shows  the  fast  convergence  of  this  algothrim.  Again,  the  Green’s 
function  for  this  case  is  uniformly  bounded  as  proven  in  Appendix  B. 

6  Conclusion 


We  have  established  the  existence  and  uniquness  of  solutions  of  the  Hertz  equation. 
A  new  numerical  method  is  devised  based  upon  the  contraction  mapping  theorem, 
and  various  examples  have  illustrated  the  usefulness  of  this  method.  For  simplicility, 
this  paper  has  dealt  solely  with  the  case  of  a.  beam.  No  clifficulities  are  encountered 
in  generalizing  the  method  for  multi-dimensional  cases  such  as  plates. 


Appendix  A 


For  the  free  vibration,  the  deflection  of  a  beam  is  goverened  by  the  PDE, 


02w  d4u> 
+ 


=  0 


0  <  x  < 


and  boundary  conditions  for  the  simply  supported  case  are 


w(  0,  t ) 
w(l,t) 


u/'(0,t)  =  0, 

=  0. 


The  corresponding  eigenvalue  problem  can  be  written  as 

W"”  =  f3W(x),  =  W"")  0  <  *  <  /, 


(AT) 


(A. 2) 


(A. 3) 


where  fi  is  an  eigenvalue  and  the  W(x)  €  L2[ 0,  /]  is  the  corresponding  eigenfunction. 
The  boundary  conditions  are 


W(0)  =  W'  (0)  =  0, 

W(l)  =  w"{l)  =  0. 

The  general  solution  can  be  expressed  as 

l,F(r)  =  c\sin(5x  +  c^cosfix  -f  c^sinliflx  +  c^coshfix, 


(A. 4) 


(A.5) 
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where  the  ci  to  C4  are  constants  to  be  determined  by  boundary  conditions  (A. 4). 
After  simplifications,  the  beam  characteristic  equation  will  be 


sin(3l  —  0  (A.G) 

The  solution  of  equation(A.G)  is  /3kl  =  &7T,  k  =  1,2, . 00 

The  normal  eigenfunctions  are 

Wk{xl)  =  \/2  sinfi^x  k  =  1,2, . 00  (A. 7) 


Lemma  A  .1  Va;,C  €  [0,/]  andVt  >  0  ,  the  Green’s  function  is  uniformly  bounded. 

Proof;  It  is  easily  checked  that  (A. 3),  (A. 4)  is  self-adjoint.  Hence,  its  Green’s 
function  can  be  expressed  as 


G{x,  C;<)  =  f^Wk(x)Wk(0SJ^H(t). 

rtt  w * 

Since , 


Wk(x) 

\Wk(x)Wk(Q\ 

|G'(;r,  (;  i)| 


< 


< 


< 


v^2 sin/3kx  V.r  6  [0,/] 

2\sin(ikxsin/3k(\  <  2 

\f:Wk(x)Wk(O^^H(t)\ 
k=  1  U’fc 


00  •  , 

Y,\mx)wk(o\\^^H(t)\ 

1 — 1 


00 


XT 


\sinwkt\ 

wk 


OO 


2H 


_1_ 


(wl  =  ,3ip/EI) 


converges,  3 M  >  0,  such  that  G(x,(j;t)  <  Af  , 


V®, C  6  [0,/]. 


(A.8) 


Appendix  B 

The  PDE  is  same  as  in  equation  (A.l),  and  the  boundary  conditions  for  the  can¬ 
tilevered  beam  are 


w(0 ,t)  =  w  (0,t)  =  0,  (B.l) 

w(l,t)  =  w  ( l,t )  =  0. 
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The  approach  to  solve  the  eigenvalue  problem  is  the  same.  Again,  we  write  the 
general  solution  and  plug  in  the  boundary  conditions  (B.l)  into  this  equation.  After 
simplification,  we  get  the  beam  characteristic  equation, 


cosftlcoshftl  =  —  1.  (B.2) 

The  orthnormal  eigenfunctions  are  determined  by  the  following  equations, 

=  j-W(x)  (B.3) 

coshftkX  —  cosftkX  sinhftkX  —  sinftkX 
coshftkl  +  cosftkl  sinhftkl  +  sinftkl 

-  f  Wl(x)dx  =  ~—J~j  A:  =  1,2, . oo. 

Jo  sm\ftkl 

Lemma  B  .1  The  orthnormal  eigenfuntions  {H4('r)}/~i  are  uniformly  bounded. 
Proof:  There  are  infinitely  many  solutions  to  the  characteristic  equation  (B.2), 
0  <  fill  <  fifil  < . <  oo,  where  fix l  =  4.73. 

Note  that  coshftkl  >  0,  coshftkX  >  0,  sinhftkl  >  0,  sinhftkX  >  0.  Vx  £  [0,/],VA:. 


Wk(x)  = 

coshftkX  -  cosftkX  sinhftkX  -  sinffox 
coshftkl  +  cos  ftkl  sinhftkl  +  sinftkl 

Ck(x)  -  Dk(x) 

(coshftkl  +  cos  ftkl  )(sinhftkl  +  sinftkl) 

where 

Ck(x)  = 

(coshftkX  -  cos ftkx)(sinh ftkl  +  sinftkl ); 

Dk(x)  = 

(sinhftkX  —  sinftkx)(coshfikl  +  cosftkl). 

Wkix) 
where  Wk{x) 
and  A'f. 


Now,  we  carry  out  the  simplifications: 


Cki'x)  —  Dk(x)  =  coshftkX  sinhftkl  —  sinhftkxcoshftkl  —  cosftkxsinhftkl  +  coshftkxsinflkl 

4-  sinftkxcoshftkl  —  sinhftkxcosftkl  +  sinfikxcosfikl  —  cosftkxsinftkl 
—  sinhftk(l  —  x)  +  sinfik(x  —  l)  —  cosftkxsinhftkl 

+  coshftkxsinftkl  +  sinftkxcoshftkl  —  sinhftkxcosftkl 


\Ck(x)-Dk(x)\  < 

< 

< 


sinhftkjl  —  x)  +  | sinftk(x  -  /)|  +  \cosfikx\sinhfikl 
+  coshftkx\sinftkl\  +  \sinftkx\coshftkl  +  sinhftkx\cosftk.l\ 
sinhftkfil  —  x)  +  1  +  sinhftkl  +  coshftkX  +  coshftkl  +  sinhftkX 
l/2e/J*(,-iC)  +  1  +  ePkl  +  ePkX  =  hk(x)  >  0; 


Since  ftkl  >  4  VAq  it  follows  that  sinhftkl  >  2,  and  coshftkl  >  2,VAr.  Hence, 


\Wl\  < 


I Ck(x)  -  Dk(x) | 

\cosIifthl  +  cosftkl\\sinhftkl  T  sinftkl \ 
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since  coshftkl  >  2; 


< 

< 

< 


_ hk(x) _ 

(coshftkl  -  1  ){sinhftkl  —  1) 

_ 2 hk(x) _ 

coshftkl(sinhftkl  —  1) 

4 hk(x) 

- - - — - —  sini 

coshftklsinhftkl 


sinhftkl  >  2; 


since  cosftklcoshftkl 

Ai 


-1  =>  cos2 fikl 

cos2ftkl  1 

sin4ftkl  Ak 


_1 _ 

cosh2 3ftkl 


and  coshftkl  >  0; 


sin2  ftklcoshftkl  =*- 


l^ik(*)|  = 


< 


< 


i  a*™ 


—  sin2  ftklcoshftkl 


4 hk{x) 


c  os  hftklsin  h  ft  k  l 


4 hk{x) 
sink  ft  kl 

8(^*(i-*)/2  +  ePkl  +  +  1) 

gPkl  —  £  — 

8{e~pkX/2  +  e^M-0  +  c~^'  +  1) 
1  —  e~2Pkl 

8(e~PlX/2  +  e^i(x-0  +  e~fhl  +  1) 

1  -  e-2/V  ’ 


Hence,  3  Mo  >  0  such  that 

|Wfc(*)|  <  Mo  V.X  e  [o,  /]:  h  =  1,2, 


For  the  cantilevered  beam,  the  differential  operator  is  also  self-adjoint.  Thus,  the 
proof  that  the  Green’s  function  is  uniformly  bounded  is  similar  to  lemma  A.l. 
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